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ABSTRACT 

For a sample of 96,951 galaxies from the Sloan Digital Sky Survey Data 
Release 3, we study the distribution of apparent axis ratios as a function of r-band 
absolute magnitude and surface brightness profile type. We use the parameter 
fracDeV to quantify the profile type (fracDeV = 1 for a pure de Vaucouleurs 
profile; f racDeV = for a pure exponential profile). When the apparent axis ratio 
9am is estimated from the moments of the light distribution, the roundest galaxies 
are very bright (M r ~ —23) de Vaucouleurs galaxies and the flattest are modestly 
bright (M r ~ —18) exponential galaxies. When the axis ratio (725 is estimated 
from the axis ratio of the 25 mag/arcsec 2 isophote, we find that de Vaucouleurs 
galaxies, at this low surface brightness, are flatter than exponential galaxies of 
the same absolute magnitude. For a given surface brightness profile type, very 
bright galaxies are rounder, on average, than fainter galaxies. We deconvolve 
the distributions of apparent axis ratios to find the distribution of the intrinsic 
short-to-long axis ratio 7, making the assumption of constant triaxiality T. For 
all profile types and luminosities, the distribution of axis ratios is inconsistent 
with a population of oblate spheroids, but is usually consistent with a population 
of prolate spheroids. Bright galaxies with a de Vaucouleurs profile (M r < —21.84, 
f racDeV > 0.9) have a distribution of g am that is consistent with triaxiality in 
the range 0.4 < T < 0.8, with mean axis ratio 0.66 < (7) < 0.69. The fainter de 
Vaucouleurs galaxies are best fit with prolate spheroids (T = 1) with mean axis 
ratio (7) ~ 0.51. 
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Subject headings: galaxies: elliptical and lenticular, cD, galaxies: fundamental 
parameters, galaxies: photometry, galaxies: spiral, galaxies: statistics 

1. INTRODUCTION 

The galaxy classification scheme of Hubble (1926) has proved durably useful. As mod- 
ified and extended by de Vaucouleurs (1959), it is still the standard method for classifying 
low-redshift galaxies with high surface brightness. The Hubble classification scheme was 
originally based on the appearance of galaxies on photographic plates. Elliptical galaxies 
have smooth elliptical isophotes; spiral galaxies have spiral arms that wind outward from a 
central bulge or bar. It was later discovered that for luminous galaxies, the surface brightness 
profile is strongly correlated with the Hubble type. If the surface brightness I is measured 
along the major axis of a galaxy's image, it is found that bright elliptical galaxies have surface 
brightness profiles that are well fit by a de Vaucouleurs, or R}l 4 law, for which log / oc — R}l 4 
(de Vaucouleurs 1948). By contrast, the azimuthally averaged surface brightness profile of a 
spiral galaxy, outside its central bulge, is typically well fit by an exponential law, log / oc —R 
(Freeman 1970). It was also eventually realized that galaxies of different Hubble type have 
different kinematic properties. The disks of spiral galaxies are rotationally flattened, with 
stars and gas on nearly circular orbits with little random motion. Bright elliptical galaxies 
{Mb < —20), by contrast, are slowly rotating and are supported mainly by their anisotropic 
velocity dispersion. 

One shortcoming of the Hubble classification scheme, imposed by necessity, is that 
elliptical galaxies are classified by their apparent two-dimensional shape, seen in projection 
on the sky, rather than their intrinsic three-dimensional shape. Consider an idealized galaxy 
whose surfaces of constant luminosity density are concentric, coaxial, similar ellipsoids, with 
principal axes of lengths a > b > c; the shape of the galaxy can then be described by 
the two axis ratios f3 = b/a and 7 = c/a. Equivalently, the shape can be described by 
the two numbers 7 and T, where the triaxiality parameter T is given by the relation T = 
(1 — /5 2 )/(l — 7 2 ). If the ellipsoidal galaxy is seen in projection, though, its isophotes will be 
concentric, coaxial, similar ellipses. The shape of the projected image can then be described 
by the single axis ratio q = B/A, where A and B are the major and minor axis length of 
any isophote. 

Although knowing the apparent axis ratio q is not, by itself, sufficient to determine 
the intrinsic axis ratios j3 and 7, the three-dimensional shape of galaxies is not beyond all 
conjecture. Two approaches to determining the three-dimensional shape of galaxies have 
been used. First, the intrinsic shape of an individual galaxy can be modeled if kinematic 
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data are available in addition to photometric data (Binney 1985; Franx, Illingworth, & de 
Zeeuw 1991). However, accurate galaxy modeling requires either multiple long-slit position 
angles (Statler 1994) or integral-field spectroscopy (Statler et al. 2004), and is best applied 
to systems with high surface brightness and large angular size. The second approach, which 
can be used in the absence of kinematic information, is to make statistical estimates of 
the distribution of f3 and 7 for a large sample of galaxies. Early estimates of the intrinsic 
shape distribution made the assumption that elliptical galaxies were oblate spheroids, with 
(3=1 (Hubble 1926; Sandage, Freeman, & Stokes 1970). More recent studies, using accurate 
surface photometry, reveal that the distribution of apparent shapes for ellipticals cannot be 
explained by a population of randomly oriented oblate spheroids; it can easily be explained, 
however, by a population of triaxial ellipsoids (Fasano & Vio 1991; Lambas, Maddox, & 
Loveday 1992; Ryden 1992; Tremblay k Merritt 1995; Alam & Ryden 2002). 

In this paper, we use the photometry-only approach to constraining the intrinsic shapes 
of galaxies, using the Sloan Digital Sky Survey Data Release 3 (SDSS DR3) as our source of 
galaxy photometry. Previous studies using the SDSS Early Data Release and Data Release 
1 studied the axis ratios of galaxies best fit by de Vaucouleurs profiles (Alam & Ryden 2002) 
and those best fit by exponential profiles (Ryden 2004). In this paper, we more closely 
examine the relation between surface brightness profile and intrinsic shape. In addition, 
we determine, for each profile type, the dependence of intrinsic shape on galaxy luminosity. 
For elliptical galaxies, the distribution of apparent shapes is known to depend on absolute 
magnitude. Elliptical galaxies with M B < —20 are rounder on average than fainter ellipticals 
(Tremblay & Merritt 1996); for a typical elliptical galaxy color of B — V rts 0.9 (Roberts & 
Haynes 1994), the dividing luminosity corresponds to My ps —20.9. In this paper, given the 
large sample size provided by the SDSS DR3, we can look for a similar shape dichotomy 
not only among galaxies with de Vaucouleurs profiles, but among galaxies with a variety of 
surface brightness profile types. 

In §2 of this paper, we describe the SDSS and the methods by which we determine the 
apparent axis ratio of the galaxies in our sample. In §3, we examine how the apparent axis 
ratios depend on the galaxy profile type and the galaxy luminosity, then use nonparametric 
kernel estimators to determine the distribution of apparent axis ratios for different samples 
of galaxies, subdivided by luminosity and profile type. In §4, we invert the distribution of 
apparent axis ratios to find the distribution of intrinsic axis ratios, assuming galaxies all have 
the same trixiality parameter T. In addition to looking at purely oblate galaxies (T = 0) 
and purely prolate galaxies (T = 1), we also examine the results for triaxial galaxies with 
T = 0.2, T = 0.4, T = 0.6, and T = 0.8. in §5, we discuss the implications of the observed 
galaxy shapes for different scenarios of galaxy formation and evolution. 
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2. DATA 

The Sloan Digital Sky Survey (York et al. 2000; Stoughton et al. 2002) will, when 
complete, provide a map of nearly one-fourth of the celestial sphere. A CCD mosaic camera 
(Gunn et al. 1998) images the sky in five photometric bands (ugriz; Fukugita et al. (1996); 
Smith et al. (2002)). The Main Galaxy Sample (MGS) of the SDSS (Stoughton et al. 2002) 
will eventually contain ~ 10 6 galaxies with r < 17.77; the mean redshift of galaxies in the 
MGS is ~ 0.1, determined from a follow-up spectroscopic survey. The SDSS Data Release 
3, issued to the astronomical community in 2004 October, contains 5282 square degrees of 
imaging data and 4188 square degrees of spectroscopic data (Abazajian et al. (2005); see 
also Stoughton et al. (2002), Abazajian et al. (2003), and Abazajian et al. (2004)). 

The SDSS DR3 data processing pipeline provides a morphological star/galaxy separa- 
tion, with extended objects being classified as 'galaxies' and point-like objects being clas- 
sified as 'stars'. For each galaxy, in each photometric band, a pair of models are fitted to 
the two-dimensional galaxy image. The first model has a de Vaucouleurs surface profile (de 
Vaucouleurs 1948): 

I(R) = J e exp (-7.67[(i?/i? e ) 1/4 - 1]) , (1) 

truncated beyond 7R e to go smoothly to zero at 8R e . The second model has an exponential 
profile: 

I(R) = I e exp (-1.68[R/R e - 1]) , (2) 

truncated beyond 3R e to go smoothly to zero at 4R e . For each model, the apparent axis 
ratio q m and the phase angle (f m are assumed to be constant with radius. The parameters 
q m , ip m , R e , and I e are varied to give the best x 2 fit to the galaxy image, after convolution 
with a double-Gaussian fit to the point spread function. 

The SDSS DR3 pipeline also takes the best-fit exponential model and the best-fit de 
Vaucouleurs model and finds the linear combination of the two that best fits the galaxy 
image. The fraction of the total flux contributed by the de Vaucouleurs component is the 
parameter fracDeV, which is constrained to lie in the interval < fracDeV < 1. The 
f racDeV parameter delineates a continuum of surface brightness profile types, from the pure 
de Vaucouleurs (fracDeV = 1) to the pure exponential (fracDeV = 0). 

The SDSS DR3 databases provide many different measures of the apparent axis ratio q 
of each galaxy in each of the five photometric bands. In this paper, we will use the r band 
data, at an effective wavelength of 6165A. (We also repeated our analysis at g and i, the 
other two bands with high signal-to-noise, and did not find significantly different results from 
our r band analysis.) Two measures of the apparent axis ratio are given by the best-fitting 
axis ratios q m for the de Vaucouleurs and exponential models. However, the algorithm for 
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fitting the models introduces quantization in the distribution of q m . Because of this artificial 
quantization, we do not use the model axis ratios as estimates of the true apparent shapes 
of galaxies. 

A useful measure of the apparent shape in the outer regions of galaxies is the axis ratio 
of the 25 mag arcsec~ 2 isophote. The SDSS DR3 data pipeline finds the best fitting ellipse 
to the 25 mag arcsec -2 isophote in each band; the semimajor axis and semiminor axis of 
this isophotal ellipse are A 25 and B 2 5- The isophotal axis ratio g 2 5 = -^25/^25 then provides 
a measure of the apparent galaxy shape at a few times the effective radius. For galaxies in 
our sample with fracDeV = 1, the mean and standard deviation of ^25/^ are 3.12 ± 1.00; 
for galaxies with fracDeV = 0, the mean and standard deviation are A2^jR e = 2.40 ± 0.36. 

Another measure of the apparent shape is g am , the axis ratio determined by the use 
of adaptive moments of the galaxy's light. The method of adaptive moments determines 
the nth order moments of a galaxy image, using an elliptical weight function whose shape 
matches that of the image (Bernstein & Jarvis 2002; Hirata & Seljak 2003). The SDSS DR3 
adaptive moments use a weight function w(x,y) that is a Gaussian matched to the size and 
ellipticity of the galaxy image I(x,y). The adaptive first order moments, 

^ = J xw(x,y)I(x,y)dxdy 
J w(x,y)I(x,y)dxdy 

and 

^ = / yw(x,y)I(x,y)dxdy 
° fw(x,y)I(x,y)dxdy 
tell us the 'center of light' of the galaxy's image. With this knowledge, we can compute the 
adaptive second order moments: 

M = J(x - x ) 2 w(x, y)I(x, y)dxdy 
fw(x,y)I(x,y)dxdy 

and so forth. The SDSS DR3 provides for each image the values of the parameters r = 
M xx + M yy , e + = (M xx — M yy )/r, and e x = 2M xy /r. The adaptive second moments can be 
converted into an axis ratio using the relation 

e\ 1/2 

' , (6) 



where e = (e+ + e 2 ) 1 / 2 . 

The adaptive moments axis ratio g am , computed in the manner given above, is not 
corrected for the effects of seeing. The SDSS DR3 also provides the fourth order adap- 
tive moments of the galaxy image, and the adaptive moments r ps f, e +iPS f, and e XiPS f of the 
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point spread function at the galaxy's location. These moments can be used to correct for 
the smearing and shearing due to seeing; such corrections are essential for studying the 
small shape changes resulting from weak lensing (Bernstein & Jarvis 2002; Hirata & Seljak 
2003). However, in this paper we will only look at well- resolved galaxies for which the seeing 
corrections are negligibly small. 

Our full sample of galaxies consists of those objects in the SDSS DR3 spectroscopic 
sample which are flagged as galaxies and which have spectroscopic redshifts z > 0.004 
(to eliminate contaminating foreground objects) and z < 0.12 (to reduce the possibility 
of weak lensing significantly distorting the observed shape). To ensure that galaxies were 
well resolved, we also required r > 7r ps f. The absolute magnitude M r of each galaxy was 
computed from the Petrosian r magnitude and the spectroscopic redshift, assuming a Hubble 
constant H = 70kms _1 Mpc^ 1 in a flat universe with f2 mi0 = 0.3 and f2 Aj0 = 0.7. No 
K-correction was applied to the data; for galaxies with normal colors at low redshift the 
K-correction in the r band is small - less than 0.2 mag for an average elliptical galaxy at 
z = 0.12, and even less for galaxies at smaller redshifts or with bluer colors (Fukugita, 
Shimasaku, & Ichikawa 1995). 

The total number of galaxies in our full sample, selected in this way, is N ga _\ = 96,951. 
Of the full sample, only 919 galaxies have fracDeV = 1, and 26,994 have fracDeV = 0. 
The remainder, constituting 71% of the full sample, are best fit by a combination of a de 
Vaucouleurs and exponential model, with < fracDeV < 1. 



3. APPARENT SHAPES 

It remains to be demonstrated that fracDeV is a useful parameter for describing the 
surface brightness profile type of galaxies. After all, a linear combination of the best de 
Vaucouleurs fit and the best exponential fit will not, in general, be the best possible de 
Vaucouleurs plus exponential fit. An alternative method of describing surface brightness 
profiles is by fitting a Sersic profile (Sersic 1968): 

log/oc-OR/it-e) 1 /", (7) 

where n can be an arbitrary number. This fitting method was used, for instance, by Blanton 
et al. (2003) in their study of SDSS galaxies. If, in fact, galaxies are well described by Sersic 
profiles, then f racDeV is a useful surrogate for the Sersic index n. Consider, for instance, 
a galaxy whose surface brightness is perfectly described by a Sersic profile of index n = 2 
and effective radius R e — X. The best-fitting exponential model for this galaxy (fitting in 
the radial region 0.1X < R < 10X) has R e = 1.10X; the best-fitting de Vaucouleurs model 
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has R e = 0.80X. Combining the models gives f racDeV = 0.51 for this n = 2 Sersic galaxy. 
A similar fit to an n = 3 Sersic galaxy yields f racDeV = 0.83. In general, if galaxies have 
Sersic profiles with 1 < n < 4, then f racDeV, as computed by the SDSS DR3 pipeline, is a 
monotonically increasing function of the Sersic index n, and thus can be used as a surrogate 
for n. 

A plot of the mean axis ratio (q) as a function of the parameter f racDeV is shown 
in Figure 1. The mean adaptive moments axis ratio (q am ), indicated by the filled circles, 
shows an increasing trend with f racDeV, from (q am ) = 0.59 for galaxies with purely expo- 
nential profiles (f racDeV = 0) to (g am ) = 0.74 for galaxies with pure de Vaucouleurs profiles 
(f racDeV = 1). However, the trend in (q) is not a linear one; for galaxies with f racDeV < 0.5, 
the value of (q) is nearly constant at (q) ~ 0.6; it is only at f racDeV > 0.5 that (q) shows an 
increasing trend with fracDeV. The mean isophotal axis ratio (#25), indicated by the open 
circles in Figure 1, shows less of a trend with f racDeV. Except in the case of nearly pure de 
Vaucouleurs profiles (fracDeV > 0.9), the axis ratio of a galaxy in its outer regions doesn't 
seem to depend on its surface brightness profile. 

For convenience in analysis, we have divided our galaxy sample into four classes, based 
on the value of fracDeV. Galaxies with fracDeV < 0.1, corresponding to a Sersic index 
n < 1.2, are called 'ex' galaxies; there are N cx = 44,289 'ex' galaxies in our sample. Galaxies 
with 0.1 < fracDeV < 0.5, corresponding to 1.2 < n < 2.0, are labeled 'ex/de' galaxies 
{N cx /dc — 36,645). Galaxies with 0.5 < fracDeV < 0.9, corresponding to 2.0 < n < 
3.3, are labeled 'de/ex' galaxies (A^e/cx — 13,780). Finally, galaxies with fracDeV > 0.9, 
corresponding to n > 3.3, are called 'de' galaxies. The small fraction of 'de' galaxies in our 
sample (iVde = 2237) is partly due to the fact that the centrally concentrated de Vaucouleurs 
galaxies are less likely to satisfy our resolution criterion, and partly due to the fact that 
galaxies with high Sersic indices are intrinsically rare; Blanton et al. (2003) estimated that 
only ~ 5% of the SDSS galaxies in their sample had Sersic index n > 3. 

The 'de' galaxies are rounder in their central regions than in their outer regions: (g am — 
Q25) = 0.083. This is consistent with the 'de' galaxies being relatively bright elliptical 
galaxies, for which the isophotal axis ratios tend to decrease with increasing semimajor axis 
length (Ryden, Forbes, & Terlevich 2001). By contrast, the 'ex' galaxies are actually slightly 
flatter in their central regions than in their outer regions: (g am — (725) = —0.017. 

For a given surface brightness profile type, there exists a relation between absolute 
magnitude and apparent axis ratio. Contour plots of mean apparent axis ratio as a function of 
both f racDeV and M r are given in Figure 2. (To give a feel for the absolute magnitude scale, 
fitting a Schechter function to the luminosity function of SDSS galaxies yields M* >r —21.4 
(Nakamura et al. 2003).) The upper panel of the figure shows that the trend in (g am ) runs 
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from the flattest galaxies at f racDeV ps and M r ps —18, where (g am ) ~ 0.52, to the roundest 
galaxies at fracDeV pa 1 and M r ps —23, where (g am ) ~ 0.83. This result is not surprising, 
since moderately bright galaxies with exponential profiles are intrinsically flattened disk 
galaxies, while extremely bright galaxies with de Vaucouleurs profiles are intrinsically nearly 
spherical giant elliptical galaxies (Tremblay & Merritt 1996). More surprising are the results 
shown in the bottom panel of Figure 2, which shows (925), the mean isophotal axis ratio. 
Here, we find that the apparently flattest galaxies, as measured by g 2 5, are not exponential 
galaxies, but galaxies with f racDeV pa 0.7 (corresponding to Sersic index n ps 2.5) and M r ps 
—20.5; these galaxies have (g 25 ) ps 0.53. The roundest galaxies, measured by (925), are not 
bright de Vaucouleurs galaxies, but bright exponential galaxies; the maximum value of (g 25 ) 
is ^ 0.74, at fracDeV ps 0, M r ps -22.5. Note also that for bright galaxies (M r < -21), the 
contours of constant (925) i n Figure 2 are nearly horizontal. That is, among bright galaxies, 
the flattening of the outer isophotes doesn't depend strongly on the surface brightness profile 
type. 

A view of the dependence of (q) on absolute magnitude for each of our four profile types, 
'ex', 'ex/de', 'de/ex', and 'de', is given in Figure 3. In the upper panel, which shows (q am ), 
note that for each profile type, there is a critical absolute magnitude M crit at which (g am ) 
is at a minimum. This critical absolute magnitude ranges from M crit ~ —20.6 for the 'de' 
galaxies to M crit ~ —19.4 for the 'ex' galaxies. At M r < M crit , the value of (g am ) increases 
relatively rapidly with increasing luminosity; at M r > M cri t, the value of (q am ) increases less 
rapidly with decreasing luminosity. At a fixed absolute magnitude, the average axis ratio of 
'de' galaxies is always greater than that of 'ex' galaxies; however, for M r < —20, the flattest 
galaxies, on average, at a given absolute magnitude are not the 'ex' galaxies, but those with 
the mixed 'de/ex' and 'ex/de' profile types. The bottom panel of Figure 3 shows (925) versus 
M r for the different profile types. In the interval —20 < M r < —22, we see that 'ex' galaxies 
have larger values of (q^) than galaxies of other profile types. 

Tremblay & Merritt (1996) divided elliptical galaxies into two classes; galaxies brighter 
than Mb ~ —20 are rounder on average than fainter ellipticals. Using a typical color for 
elliptical galaxies of B — V ~ 0.9, this corresponds to an absolute magnitude in the r band 
of M r . Pi -21.2, using the transformation r = B- 1.44(5 — V) + 0.12 (Smith et al. 2002). In 
our data, we do not see an abrupt jump in (g am ) for 'de' galaxies at M r = —21.2, but rather 
a more gradual increase between M r ps —20.6 and M r Pi —22.7. In addition, it is true for 
any profile type that highly luminous galaxies are rounder, on average, than less luminous 
galaxies. 

Knowing the mean shape as a function of M r for each profile type ('de', 'de/ex', 'ex/de', 
'ex') does not tell us the complete distribution of axis ratios, f{q). In order to compare the 
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distribution of axis ratios for highly luminous galaxies with that for less luminous galaxies, 
we determined, for each profile type, the dividing absolute magnitude M such that the 
distribution f(q am ) for galaxies with M r < M is maximally different from the distribution 
/(<Zam) f° r galaxies with M r > M . We measure the difference between the distributions 
by applying a Kolmogorov-Smirnov (K-S) test, and finding the value of Mq that minimizes 
the K-S probability Prs- The values of the dividing magnitude M for each profile type, 
and the associated probability P KS , are given in Table 1. (Using a Student's t-test for the 
significance of the difference of the means yields similar values of M , differing by less than 
0.2 mag from the values found using the K-S test.) 

To estimate the distribution of apparent axis ratios f(q) for each of our galaxy subsam- 
ples, we use a nonparametric kernel technique (Vio et al. 1994; Tremblay & Merritt 1995; 
Ryden 1996). Given a sample of N axis ratios, q±, q2, ■ ■ ., qN, the kernel estimate of the 
frequency distribution f(q) is 



N 

Q-Qi 



(8) 



hi 

where K(x) is the kernel function, normalized so that 

K(x)dx = 1 . (9) 



To ensure that our estimate / is smooth, we use a Gaussian kernel, with K(x) oc exp(— x 2 /2). 
The parameter hi in equation (8) is the kernel width; too small a width introduces noise into 
the estimate, while too great a kernel width produces excessive smoothing. We choose the 
kernel width hi by using the standard adaptive two-stage estimator of Abramson (1982). In 
this technique, an initial estimate f is made using a global fixed kernel width 

h = 0.9A/N - 2 , (10) 

with A = min(<T, Q4/I.34), where is the standard deviation and Q4 is the interquartile 
range of the axis ratios. For samples that are not strongly skewed, this formula minimizes 
the mean square error (Silverman 1986; Vio et al. 1994). The final estimate / is then 
computed, using at each data point qi the kernel width 



hi = h 



G/h{qi) V2 , (11) 



where G is the geometric mean of f over all values of i. 

The axis ratio for a galaxy must lie in the range < q < 1. To ensure that / = for 
q < and q > 1, we impose reflective boundary conditions at q = and q — 1. In practice, 
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this means replacing the simple Gaussian kernel with the kernel 

'-'(x)^)^^)' (12) 

Although use of reflective boundary conditions ensures that / remains within bounds, it 
imposes the possibly unphysical constraint that df/dq = at q = and q — 1. 

The estimated distributions /(gam) are shown in the left column of Figure 4, for the 
'de', 'de/ex', 'ex/de', and 'ex' galaxies; the equivalent distributions /(?2s) f° r the isophotal 
shapes are shown in the right column. In each panel of Figure 4, the heavier line indicates 
the shape distribution for the bright galaxies, and the lighter line indicates the distribution 
for the fainter galaxies. The dotted lines indicate the 98% error intervals found by bootstrap 
resampling. In our bootstrap analysis, we randomly selected iV data points, with substi- 
tution, from the original sample of N axis ratios. We then created a new estimate / from 
the bootstrapped data. After doing 1000 bootstrap estimates, we then determined the 98% 
error intervals: in Figure 4, 1% of the bootstrap estimates lie above the upper dotted line, 
and 1% of the bootstrap estimates lie below the lower dotted line. 

Note from the upper left panel of Figure 4 that the bright 'de' galaxies have a distribution 
that peaks at g am 0.84, similar to the result found for relatively bright elliptical galaxies 
(Benacchio & Galletta 1980; Fasano & Vio 1991; Franx, Illingworth, & de Zeeuw 1991; 
Lambas, Maddox, & Loveday 1992). The faint 'de' galaxies have a flatter modal shape 
(<?am ~ 0.70) and a broader distribution of shapes. Both the bright and faint 'de' galaxies 
have a relative scarcity of galaxies with g am = 1; this is the usual sign that the galaxies 
in the population cannot be purely oblate. For 'de' galaxies, the difference between the 
adaptive moments shapes and the isophotal shapes is intriguing. Note that /(^s) for the 
'de' galaxies, shown in the upper right panel of Figure 4, peaks at (725 ~ 0.76 for the bright 
galaxies. For the fainter 'de' galaxies, the distribution of isophotal axis ratios, though it 
peaks at g 25 ~ 0.74, is very broad. An axis ratio distribution / that is nearly constant 
over a broad range of q is the signature of a population of flattened disks. The distribution 
/(g 25 ) drops at q < 0.2, indicating that the disks are not infinitesimally thin, and at q > 0.9, 
indicating that the disks are not perfectly circular. 

The distribution of shapes for 'de/ex' galaxies are seen in the second row of Figure 4. 
As with all profile types, the bright galaxies are rounder on average than the fainter galaxies. 
Both bright and faint 'de/ex' galaxies have a scarcity of nearly circular galaxies. Particularly 
noteworthy is the distinct difference between f(q am ) for the bright 'de' galaxies, very strongly 
peaked at g am ~ 0.84, and the equivalent distribution /(g am ) for the bright 'de/ex' galaxies, 
which is much broader and peaks at the flatter shape of q am ~ 0.67. 



- 11 - 



The 'ex/de' and 'ex' galaxies in our sample have, by definition, fracDeV < 0.5 and 
hence Sersic index n < 2. We expect these galaxies to be predominantly disk-dominated 
spiral galaxies. Determining the true distribution of intrinsic axis ratios for spiral galaxies 
is complicated by the presence of dust. A magnitude-limited sample of dusty spiral galaxies 
will show a deficit of galaxies with high inclination and low apparent axis ratio. This deficit 
is due to internal extinction by dust; in the B band, for instance, a spiral galaxy of type Sc 
will be up to 1.5 mag fainter when seen edge-on than when seen face-on (Huizinga & van 
Albada 1992). In the r band, the inclination-dependent dimming is not as great, but will 
still affect the observed distribution of axis ratios in a magnitude-limited sample such as the 
SDSS spectroscopic galaxy sample (Ryden 2004). A dusty disk galaxy which would qualify 
as a bright 'ex' galaxy if seen face-on might only qualify as a faint 'ex' galaxy when seen 
edge-on. Similarly, a less luminous galaxy, which would qualify only as a faint 'ex' galaxy 
when seen face-on, might drop out of our sample entirely when seen edge-on. The main 
effect of dust on the distribution of q for a magnitude-limited survey of spiral galaxies is 
to reduce the number of galaxies with q < 0.5; the distribution for q > 0.5 is not strongly 
affected (Huizinga & van Albada 1992). 

With the caveat that we are undercounting edge-on, low-g galaxies, the observed axis 
ratios for 'ex/de' and 'ex' galaxies are shown in the bottom two rows of Figure 4. The shapes 
of 'ex/de' galaxies, as measured by either q am or g 2 5, spread over a wide range of axis ratios, 
with bright 'ex/de' galaxies being slightly rounder, on average, than the fainter galaxies. 
The shapes of 'ex' galaxies are very similar to the shapes of 'ex/de' galaxies. For the 'ex/de' 
and 'ex' galaxies, the dip in / at q > 0.9 is a consequence of the fact that the disks of these 
galaxies are not perfectly circular. 



4. INTRINSIC SHAPES 

Having a smooth kernel estimate / for the distribution of apparent axis ratios q permits 
us to estimate the distribution of intrinsic axis ratios 7, given the assumption that all galaxies 
have the same value of the triaxiality parameter T. We don't necessarily expect all the 
galaxies of a given profile type to have the same triaxiality; cosmological simulations, as well 
as simulations of merging galaxies, generally give a range of T (see Statler et al. (2004) and 
references therein) for bright galaxies. However, by assuming uniform triaxiality and seeing 
which values of T yield physically plausible results, we can get a feel for which values of T 
are data- friendly, and which are incompatible with the data. Let At (7) be the estimated 
distribution of intrinsic axis ratios, given an assumed value for T. The relation between f(q) 
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and N T (~f) is (Binney 1985): 

/(g) = r P T (qh)N T (^ , (13) 

where -Pr(g|7)<^g is the probability that a galaxy with intrinsic short-to-long axis ratio 7 and 
triaxiality T has an apparent axis ratio in the range q — > q + dq, when viewed from a random 
orientation. 

Equation (13) represents a Volterra equation of the first kind which can be solved by 
standard numerical techniques. We first compute / on 200 points equally spaced in q, 
using the kernel technique previously described. We then compute Pt on an 200 x 200 
grid in (g, 7) space. For each value of 7, we compute an adaptive kernel estimate for Pt 
by randomly choosing a large number of viewing angles (~ 10 5 ), then computing what the 
apparent axis ratio q would be for each viewing angle, given the assumed values of 7 and T 
(Binney 1985). The kernel estimate of the distribution of q is then our estimate of Pr(g|7)- 
Equation (13), in its discretized form, can be written 

k 

//.. Yl /, '-^ A >-/ • ( 14 ) 

Since Pr,kj = for j > k (the apparent axis ratio of a galaxy can't be less than its intrinsic 
axis ratio), the matrix P T is lower triangular, and equation (14) can be inverted simply by 
forward substitution (Press et al. 1992). 

We performed the inversion to find Nt using six different values of the trixiality: T = 
(oblate), T = 0.2, T = 0.4, T = 0.6, T = 0.8, and T = 1 (prolate). When our ultimate 
goal is finding N T , determining the optimal kernel width h for / is a vexed question. It is 
generally true that the best value of h for estimating Nt is greater than the best value of 
h for purpose of estimating / (Tremblay & Merritt 1995, 1996); Nt is a deconvolution of 
/, and the deconvolution process increases the noise present. Tremblay & Merritt (1995) 
proposed, for the purposes of estimating Nt, using the value of h that minimized 

BMISE(fc) = / [JV T (7) - N*(-f)] 2 d-f , (15) 
Jo 

where N T is an estimate of N T from a bootstrap resampling of the original data. Un- 
fortunately, we found, as Tremblay & Merritt (1995) did, that BMISE(/i) returns badly 
oversmoothed estimates of N T - Thus, we fell back on using our subjective impressions of the 
smoothness of N T , and ended by taking the value of h given by equation (10) and multiplying 
by a factor of 1.5. 
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To perform the inversion and find the estimated distribution of 7, we assumed that all 
galaxies have the same triaxiality T. If this assumption is incorrect, then the inversion of 
equation (14) may result in < for some range of 7, which is physically impossible. 
To exclude our hypothesis (that all galaxies have triaxiality T) at some fixed statistical 
confidence level, we can repeat the inversion for a large number of bootstrap resamplings 
of the original data set. In practice, we did 1000 resamplings, and used them to place 98% 
confidence intervals on N T . If the confidence interval falls below zero for some valus of 7, 
then the hypothesized value of T can be rejected at the 99% (one-sided) confidence level. 

Figure 5 shows the distribution of intrinsic shapes for 'de' galaxies (fracDeV > 0.9), 
when the adaptive moments axis ratio q am is used as the estimator of the apparent shape. 
The heavy line is the distribution for the bright 'de' galaxies (M r < —21.8) and the light 
line is the distribution for the fainter 'de' galaxies. The assumed value of T is shown in each 
panel, ranging from purely oblate shapes, with T = 0, in the upper left corner to purely 
prolate shapes, with T = 1, in the lower right corner. 

An eyecatching feature of Figure 5 is the oscillatory nature of N T in the triaxial cases. 
The multiple peaks in N T at high 7 result from the shape of the conditional probability 
function Pt(?|t) for highly triaxial galaxies. When galaxies are axisymmetric (T = or 1), 
the conditional probability peaks at q — 7, the minimum possible apparent axis ratio. For 
triaxial galaxies, however, the conditional probability has one or more local maxima at q > 7. 
To see why this can result in oscillatory solutions consider for faint 'de' galaxies when 
T = 0.8; this is shown by the light line in the lower left panel of Figure 5. The distribution 
has a local maximum at 7 = 0.42. For T = 0.8, the conditional probability Pt(?|7 — 0.42) 
peaks at q — 0.57. Thus, the large number of 7 « 0.42 galaxies result in many galaxies with 
an apparent axis ratio q ~ 0.57; so many, in fact, that N T must be made very small in the 
range 0.42 < 7 < 0.57 in order to avoid overproducing q « 0.57 galaxies. In fact, as seen in 
Figure 5, N T has a local minimum at 7 = 0.54, at which N T falls slightly below zero. Since 
Pr(<7 17 = 0.54) peaks at q — 0.65, this produces a deficit of q ~ 0.65 galaxies, which must 
be compensated for by making very large in the range 0.54 < q < 0.65. In fact, Nt 
has another local maximum at q — 0.62. And so the oscillations continue, with decreasing 
wavelength, until 7 = 1 is reached. Tremblay & Merritt (1995) found that ^(7) for bright 
elliptical galaxies was significantly bimodal if the galaxies were assumed to be highly triaxial; 
this bimodality also had its origin in the shape of the conditional probability function for 
triaxial galaxies. 

For both bright and faint 'de' galaxies, using q am as the apparent axis ratio, oblate fits 
are statistically unacceptable, producing a negative number of galaxies with 7 > 0.9 (see the 
upper left panel of Figure 5). For bright 'de' galaxies, statistically acceptable fits are found 
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for T = 0.4, 0.6, and 0.8; that is, their 98% confidence intervals never fall entirely below 
zero. The mean intrinsic axis ratio for the acceptable fits ranges from (7) = 0.66 for T = 0.4 
to (7) = 0.69 for T = 0.8. The permissible fits for the bright 'de' galaxies are consistent with 
the deduced triaxial shape of the nearby bright (M r ~ —22.2) elliptical galaxy NGC 4365, for 
which a combination of photometric and kinematic data yields T ~ 0.45 and 7 ~ 0.6 (Statler 
et al. 2004). Although the bright 'de' galaxies are best fit by highly triaxial shapes, the faint 
'de' galaxies are best fit by nearly prolate shapes, with T = 0.8 and T = 1 both giving 
statistically acceptable fits. (The highly oscillatory solution for T = 0.8 may be physically 
dubious - why should galaxy shapes be quantized? - but it is statistically acceptable.) The 
mean intrinsic shape of faint 'de' galaxies is (7) = 0.53 if T = 0.8 and (7) = 0.58 if T — 1. 

Figure 6 shows the deduced distribution of intrinsic shapes for 'de' galaxies when the 
isophotal axis ratio g 2 5 is used, rather than the adaptive moments axis ratio. The bright 'de' 
galaxies have a statistically acceptable fit when T = 0.4, T = 0.6, T = 0.8, and T = 1. The 
mean intrinsic shape ranges from (7) = 0.51 when T = 0.4 to (7) = 0.62 when T = 1. The 
fainter 'de' galaxies are acceptably fit, assuming constant triaxiality, only when T = 1, which 
results in (7) = 0.51. The T = case, which can be rejected at the 99% confidence level but 
not at the 99.9% level, would produce (7) = 0.28; this axis ratio is flatter than that of an 
ice hockey puck. (Although the T = and T = 0.2 fits are statistically unacceptable at the 
99% confidence level, the data are consistent with a population of nearly oblate shapes if we 
relax the assumption of uniform T. For instance, the faint 'de' galaxies can be fit at the 99% 
confidence level, with a Gaussian distribution of T peaking at T = and with a T = 0.2.) 

The shapes of 'de/ex' galaxies can be analyzed in the same way as the shapes of 'de' 
galaxies. For instance, the distributions of intrinsic shapes for 'de/ex' galaxies, using g am 
as the apparent axis ratio, are shown in Figure 7. For the brighter galaxies with 'de/ex' 
profiles, acceptable fits are found when T = 0.8, yielding (7) = 0.51, and when T = 1, 
yielding (7) = 0.57. For the fainter 'de/ex' galaxies, only the T = 1 fit is statistically 
acceptable, yielding a mean intrinsic axis ratio of (7) = 0.48 for the prolate galaxies. Notice 
that bright 'de/ex' galaxies, faint 'de/ex' galaxies, and faint 'de' galaxies all are consistent 
with having shapes, at least in their inner regions, which are prolate, or nearly prolate, with 
a typical axis ratio of (7) ~ 1/2. Bright 'de' galaxies are distinctly different, with shapes 
that are consistent with their being highly triaxial, with a typical axis ratio of (7) ~ 2/3. 

The intrinsic shape distributions for 'de/ex' galaxies, using g 2 5 as the apparent shape 
measure, are shown in Figure 8. For both the bright and faint 'de/ex' galaxies, the T = 1 
fit is statistically acceptable. For the prolate fits, the average shape of the bright 'de/ex' 
galaxies is (7) = 0.55, and the average shape of the faint 'de/ex' galaxies is (7) = 0.47. Thus, 
given the prolate hypothesis, 'de/ex' galaxies are only slightly flatter in their outer regions, 
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whose shape is measured by #25, than in their inner regions, whose shape is measured by 

With the repeated caveat that we are undercounting edge-on, \ow-q galaxies, the results 
for the 'ex/de' galaxies are presented in Figures 9 and 10. When q am is used as the measure 
of apparent shape (Figure 9), the bright 'ex/de' galaxies have a statistically acceptable fit 
when T = 0.8, which results in (7) = 0.49; the T = 1 fit, which goes slightly negative at 
7 > 0.94, yields a mean intrinsic axis ratio of (7) = 0.55. None of the six tested triaxialities 
gave an acceptable fit, at the 99% confidence level, to the faint 'ex/de' galaxies. The best (or 
the 'least bad') of the fits, with T = 1, yields (7) = 0.49. When ^25 is used as the measure 
of apparent shape (Figure 10), none of the tested values of T gives a fit acceptable at the 
99% confidence level. The best fit, for both bright and faint galaxies, is given by T = 1, for 
which (7) = 0.55 for the bright 'ex/de' galaxies and (7) = 0.50 for the faint 'ex/de' galaxies. 

As with the 'ex/de' galaxies, the sample of 'ex' galaxies is affected by an undercount 
of low-g galaxies. With this caveat, the results for the 'ex' galaxies are shown in Figures 11 
and 12. In general, the results are similar to those for 'ex/de' galaxies. When the axis ratios 
of 'ex' galaxies are estimated using g am (Figure 11), the bright 'ex' galaxies are acceptably 
fit when T = 0.8, yielding (7) = 0.49; the T = 1 fit for the bright 'ex' galaxies, which goes 
slightly negative at high 7, yields (7) = 0.56. The best fit to the faint 'ex' galaxies is provided 
by the T = 1 inversion, which yields (7) = 0.46. When ^25 is used as the estimate of the 
apparent axis ratio (Figure 12), the results are extremely similar. For bright 'ex' galaxies, 
only T = 0.8, of the tested triaxialities, gave a fit acceptable at the 99% confidence level, 
yielding (7) = 0.51; the not-as-good T = 1 fit gave (7) = 0.57. The best fit for faint 'ex' 
galaxies had T = 1, (7) = 0.49. 

Although the best fits to the 'ex/de' and 'ex' galaxies, under the assumption of con- 
stant triaxiality, were prolate or nearly prolate, we have an a priori knowledge that bright 
galaxies with Sersic index n < 2 are generally rotationally supported disks. The SDSS DR3 
photometry of 'ex/de' and 'ex' galaxies is consistent with their being nearly oblate objects 
only if they have a range of T, instead of being forced into a straitjacket of uniform T. 

5. DISCUSSION 

Constraining the flattening 7 and triaxiality T of galaxies with different luminosity and 
profile type provides potentially useful clues for the study of galaxy formation and evolution. 
For example, the dissipationless merger of two equal-mass disk galaxies embedded within 
dark halos produces a merger remnant whose luminous component has an approximate de 



-16- 



Vaucouleurs profile (Barnes 1992). The values of 7 and T of the remnant depend on the 
initial orbital parameters of the merging galaxies, but the remnant is typically quite flat 
and prolate. For a set of eight simulations with different initial orbits, Barnes (1992) found 
(7) ~ 0.5 and (T) ~ 0.7 for the most tightly bound quartile of the luminous particles (the 
inner region of the galaxy). This is significantly flatter and more prolate than the shape we 
deduced from the distribution of g am for bright 'de' galaxies. The oblateness of the merger 
remnant is increased, however, if the initial merging galaxies are unequal in mass (Barnes 
1998; Naab & Burkert 2003). 

Including gaseous dissipation into a simulated merger can strongly affect 7 and T. 
As the gas dissipates and falls toward the center of the merger remnant, the central mass 
concentration will scatter stars on low angular momentum box orbits, thus acting to increase 
7 and decrease T. For instance, a simulated remnant that has 7 p=* 0.55 and T ps 0.8 in the 
absence of dissipation will have 7 ps 0.65 and T ps 0.5 if gas dissipation is added (Barnes 
1998). The deduced shapes of bright 'de' galaxies are consistent with their being merger 
remnants, as long as the mergers were dissipational. 

In the standard hierarchical clustering model for structure formation, however, we do 
not expect most elliptical galaxies to have formed by the relatively recent merger of a pair of 
spiral galaxies. More typically, an elliptical galaxy will have formed by the successive merger 
of a number of smaller galaxies or subgalactic "clumps" . Simulations of multiple mergers in 
small groups reveal that the final merger remnant generally has a surface brightness profile 
similar to a de Vaucouleurs law, but that the shape of the merger remnant depends on the 
assumed initial conditions (Weil & Hernquist 1996; Garijo, Athanassoula, & Garcia-Gomez 
1997; Zhang et al. 2002). For instance, Zhang et al. (2002) simulated dissipationless mergers 
of dense clumps embedded in a smooth dark halo. The profile of the merger remnant, in 
all cases, was well fit by a de Vaucouleurs law. The triaxiality of the remnant depended 
primarily on whether the system was initially in a state of virial equilibrium or in a state of 
non-virial collapse. The virial case produced more oblate remnants ((T) ~ 0.3), while the 
collapse case produced more prolate remnants ((T) ~ 0.7). The virial case also produced 
less flattened remnants: (7) ~ 0.8 for the merger remnants of initially virialized systems 
versus (7) ~ 0.6 for initially collapsing systems. The value of 7 also depended on the initial 
clump-to-halo mass ratio; if the clumps contributed only a small fraction of the total mass, 
they produced a flatter remnant than if they contributed all the mass, with no dark halo. 
The best fitting parameters for the bright 'de' galaxies, /i 7 = 0.66 and T = 0.43, suggest 
that they are too flattened and too triaxial to have been formed primarily by mergers within 
a virialized group, which produce (7) ~ 0.8 and (T) ~ 0.3. (Remember, the effects of 
dissipation will only act to increase 7 and decrease T.) 
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Another influence that is capable of increasing 7 and decreasing T is the presence of a 
central supermassive black hole. A central mass concentration affects the structure of the 
surrounding galaxy by disrupting box orbits (Gerhard & Binney 1985; Norman, May, & van 
Albada 1985). The resulting chaotic orbits lead to an altered shape for the galaxy; the trend 
is toward shapes that are nearly spherical in the inner regions and nearly oblate in the outer 
regions (Merritt & Quinlan 1998). In the inner regions, the orbits can support a shape that 
is nearly oblate, T ~ 0.25, or strongly triaxial, T ~ 0.5, but not nearly prolate, T ~ 0.75 
(Poon & Merritt 2004). The time scale for the shape evolution depends on the mass of the 
central black hole relative to the total stellar mass of the galaxy. When the central black 
hole has a mass equal to 1% of the mass in stars, its shape at the half-mass radius evolves 
in ~ 40 times the half-mass orbital period; when the mass is only 0.3% of the stellar mass, 
the evolution time scale increases to ~ 200 times the orbital period. The observed relation 
between central black hole mass and bulge mass (Merritt & Ferrarese 2001; Marconi & Hunt 
2003; Haring & Rix 2004) predicts that an elliptical galaxy will typically have a central black 
hole with a mass < 0.2% of the total mass in stars. For such relatively low black hole masses, 
the shape evolution time scale at the half-mass radius and beyond will exceed the age of the 
galaxy. Thus, a central black hole will not generally affect the isophotal axis ratio g 25 . 

The isophotal shape, g 25 , measures the apparent shape of a galaxy at a few times 
the effective radius. Converted to physical units, the average semimajor axis A 2 $ of the 
isophotal ellipse ranges from (A 25 ) = 2.4 kpc at M r = —16 to (A 25 ) = 45kpc at M r = 
—23. At relatively high luminosities (M r < —20.5), the mean value of g 25 increases with 
luminosity; bright galaxies are rounder than fainter galaxies at their outer regions. We 
have also uncovered the puzzling result that at a given luminosity, the mean value of g 2 5 is 
larger for 'ex' galaxies than for 'de' galaxies. Why should 'ex' galaxies, which are flattened 
disks in their central regions, be rounder than 'de' galaxies in their outskirts? At such great 
distances from the galaxy's center, we expect that tidal distortions from neighboring galaxies 
will significantly affect the isophotal shape. Galaxies in dense regions, such as the cores 
of clusters, are subject to more interactions with neighboring galaxies, and may well differ 
systematically in shape from field galaxies, which are less frequently harassed. Since elliptical 
galaxies are found preferentially in dense environments, this conjectured environment - shape 
relation would then translate into a relation between profile type and shape. A future paper, 
studying the relation between apparent axis ratio and environment for SDSS DR3 galaxies 
(Kuehn & Ryden, 2005, in preparation), will test this conjecture. 

Richard Pogge, Fred Kuehn, and the anonymous referee provided useful assistance and 
comments. 

Funding for the creation and distribution of the SDSS Archive has been provided by 



-18- 



the Alfred P. Sloan Foundation, the Participating Institutions, the National Aeronautics 
and Space Administration, the National Science Foundation, the U.S. Department of En- 
ergy, the Japanese Monbukagakusho, and the Max Planck Society. The SDSS website is 
http://www.sdss.org/. The SDSS is managed by the Astrophysical Research Consortium 
(ARC) for the Participating Institutions. The Participating Institutions are The University 
of Chicago, Fermilab, the Institute for Advanced Study, the Japan Participation Group, The 
Johns Hopkins University, Los Alamos National Laboratory, the Max-Planck-Institute for 
Astronomy (MPIA), the Max-Planck- Institute for Astrophysics (MPA), New Mexico State 
University, University of Pittsburgh, Princeton University, the United States Naval Obser- 
vatory, and the University of Washington. 

REFERENCES 

Abazajian, K., et al. 2003, AJ, 126, 2081 
Abazajian, K. et al. 2004, AJ, 128, 502 

Abazajian, K. et al. 2005, AJ, submitted (astro-ph/0410239) 
Abramson, I. S. 1982, Annals of Statistics, 10, 1217 
Alam, S. M. K., & Ryden, B. S. 2002, ApJ, 570, 610 
Barnes, J. E. 1992, ApJ, 393, 484 

Barnes, J. E. 1998, in Galaxies: Interactions and Induced Star Formation, ed. D. Friedli, L. 
Martinet, & D. Pfenniger (Saas-Fee Adv. Course 26; Berlin: Springer), 275 

Benacchio, L., & Galletta, G. 1980, MNRAS, 193, 885 

Bernstein, G. M., & Jarvis, M. 2002, A J, 123, 583 

Binney, J. 1985, MNRAS, 212, 767 

Blanton, M. R., et al. 2003, ApJ, 594, 186 

de Vaucouleurs, G. 1948, Ann. d'Astrophys. 11, 247 

de Vaucouleurs, G. 1959, Handb. Phys. 53, 275 

Fall, M., & Frenk, C. S. 1983, A J, 88, 1626 

Fasano, G., & Vio, R. 1991, MNRAS, 249, 629 



-19- 

Franx, M., Illingworth, G., & de Zeeuw, T. 1991, ApJ, 383, 112 
Freeman, K. C. 1970, ApJ, 160, 811 

Fukugita, M., Shimasaku, K., Ichikawa, T. 1995, PASP, 107, 945 

Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shimasaku, K., & Schneider, D. P. 1996, 
AJ, 111, 1748 

Garijo, A., Athanassoula, E., & Garcia-Gomez, C. 1997, A&A, 327, 930 

Gerhard, O. E., & Binney, J. 1985, MNRAS, 216, 467 

Gunn, J. E., et al. 1998, AJ, 116, 3040 

Haring, N., & Rix, H.-W. 2004, ApJ, 604, L89 

Hirata, C, & Seljak, U. 2003, MNRAS, 343, 459 

Hubble, E. 1926, ApJ, 64, 321 

Huizinga, J. E., & van Albada, T. S. 1992, MNRAS, 254, 677 
Lambas, D. G., Maddox, S., & Loveday J. 1992, MNRAS, 258, 404 
Marconi, A., & Hunt, L. K. 2003, ApJ, 589, L21 
Merritt, D., & Ferrarese, L. 2001, MNRAS, 320, L30 
Merritt, D., & Quinlan, G. D. 1998, ApJ, 498, 625 
Naab, T., & Burkert, A. 2003, ApJ, 597, 893 

Nakamura, O., Fukugita, M., Yasuda, N., Loveday, J., Brinkmann, J., Schneider, D. P., 
Shimasaku, K., SubbaRao, M. 2003, AJ, 125, 1682 

Norman, C. A., May, A., & van Albada, T. S. 1985, ApJ, 296, 20 

Poon, M. Y., & Merritt, D. 2004, ApJ, 606, 774 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical Recipes: 
The Art of Scientific Computing, Second Edition (Cambridge: Cambridge University 
Press), ch. 2 

Roberts, M. S., & Haynes, M. P. 1994, ARA&A, 32, 115 



-20- 

Ryden, B. S. 1992, ApJ, 396, 445 
Ryden, B. S. 1996, ApJ, 461, 146 
Ryden, B. S. 2004, ApJ, 601, 214 

Ryden, B. S., Forbes, D. A., & Terlevich, A. I. 2001, MNRAS, 326, 1141 

Sandage, A., Freeman, K. C, & Stokes, N R. 1970, ApJ, 160, 831 

Sersic, J. L. 1968, Atlas de Galaxias Australes (Cordoba: Obs. Astron.) 

Silverman, B. W. 1986, Density Estimation for Statistics and Data Analysis (New York: 
Chapman & Hall) 

Smith, J. A., et al. 2002, AJ, 123, 2121 

Statler, T. S. 1994, ApJ, 425, 458 

Statler, T. S., Emsellem, E., Peletier, R. F., & Bacon, R. 2004, MNRAS, 353, 1 

Stoughton, C, et al. 2002, AJ, 123, 485 

Tremblay, B., & Merritt, D. 1995, AJ, 110, 1039 

Tremblay, B., & Merritt, D. 1996, AJ, 111, 2243 

Weil, M. L., & Hernquist, L. 1996, ApJ, 460, 101 

Vio, R., Fasano, C, Lazzarin, M., & Lessi, O. 1994, A&A, 289, 640 

York, D. G. et al. 2000, AJ, 120, 1579 

Zhang, B., Wyse, R. F. C, Stiavelli, M., & Silk, J. 2002, MNRAS, 332, 647 



This preprint was prepared with the AAS IATgX macros v5.2. 



- 21 - 



Table 1. Dividing Line Between Bright and Faint Galaxies 



Profile Type 


M 


N(M r < M ) 


N(M r > M ) 


log P K s 


de 


-21.84 


926 


1311 


-52.9 


de/ex 


-21.35 


5235 


8545 


-109.2 


ex/de 


-21.00 


11340 


25305 


-118.6 


ex 


-20.48 


17101 


27188 


-301.9 




.2 .4 .6 .8 1 

fracDeV 

Fig. 1. — Mean apparent axis ratio (q) as a function of the SDSS fitting parameter f racDeV. 
Filled circles indicate the axis ratio estimated using the adaptive moments technique; open 
circles indicate the axis ratio of the 25 mag/arcsec 2 isophote. The error bars on the filled 
circles indicate the estimated error of the mean; the error bars on the open circles (not shown 
for clarity) are of similar size. 
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fracDeV 

Fig. 2. — Upper panel: Mean apparent axis ratio g am , estimated using adaptive moments, as a 
function of absolute magnitude M r and fitting parameter f racDeV. Averages were computed 
in bins of width 0.25 mag in M r and 0.05 in fracDeV; the heavy line excludes the region 
with fewer than 25 galaxies per bin. Contours are drawn at q — 0.6, 0.7, and 0.8 (solid lines) 
and at q = 0.55, 0.65, and 0.75 (dotted lines). Lower panel: Same as the upper panel, using 
the isophotal axis ratio #25 as the shape estimate. 
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Fig. 3. — Upper panel: Mean apparent axis ratio g am , estimated using the adaptive moments 
technique, as a function of absolute magnitude. The four different galaxy profile types (de, 
de/ex, ex/de, ex) are described in the text. Lower panel: Same as the upper panel, using 
the apparent axis ratio g 2 5 of the 25 mag/arcsec 2 isophote. 
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Fig. 4. — Le/t column: Distribution of the adaptive moments axis ratio g am . From top to 
bottom, the results are shown for profile types 'de', 'de/ex', 'ex/de', and 'ex'. In each panel, 
the heavy line is the distribution for bright galaxies, and the light line is the distribution 
for faint galaxies. The dotted lines indicate the 98% confidence interval found by bootstrap 
resampling. The horizontal error bars indicate the kernel width h. Right column: The same 
as the left column, but using the 25 mag/arcsec 2 isophotal axis ratio, g 25 , as the shape 
estimator. 
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Fig. 5. — Distribution of intrinsic axis ratios for galaxies of profile type 'de' (f racDeV > 0.9), 
using the adaptive moments estimate of q. The heavy line is the distribution for galaxies 
with M r < —21.8; the light line is for galaxies with M r > —21.8. The dotted lines indicate 
the 98% confidence intervals, estimated by bootstrap resampling. The assumed value of the 
triaxiality T is given in each panel. 
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Fig. 6. — Same as Figure 5, but using the apparent axis ratio of the 25 mag/arcsec 2 isophote 
as the shape estimate. 




Fig. 7. — Distribution of intrinsic axis ratios for galaxies of profile type 'de/ex' (0.5 < 
f racDeV < 0.9), using the adaptive moments estimate of q. The heavy line is the distribution 
for galaxies with M r < —21.35; the light line is for galaxies with M r > —21.35. The dotted 
lines indicate the 98% confidence intervals, estimated by bootstrap resampling. The assumed 
value of the triaxiality T is given in each panel. 




Fig. 8. — Same as Figure 7, but using the apparent axis ratio of the 25 mag/arcsec 2 isophote 
as the shape estimate. 
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Fig. 9. — Distribution of intrinsic axis ratios for galaxies of profile type 'ex/de' (0.1 < 
f racDeV < 0.5), using the adaptive moments estimate of q. The heavy line is the distribution 
for galaxies with M r < —21.0; the light line is for galaxies with M r > —21.0. The dotted 
lines indicate the 98% confidence intervals, estimated by bootstrap resampling. The assumed 
value of the triaxiality T is given in each panel. 
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Fig. 10. — Same as Figure 9, but using the apparent axis ratio of the 25 mag/arcsec 2 isophote 
as the shape estimate. 
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Fig. 11. — Distribution of intrinsic axis ratios for galaxies of profile type 'ex' (f racDeV < 0.1), 
using the adaptive moments estimate of q. The heavy line is the distribution for galaxies 
with M r < —20.5; the light line is for galaxies with M r > —20.5. The dotted lines indicate 
the 98% confidence intervals, estimated by bootstrap resampling. The assumed value of the 
triaxiality T is given in each panel. 
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Fig. 12. — Same as Figure 11, but using the apparent axis ratio of the 25 mag/arcsec 2 
isophote as the shape estimate. 



